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ABSTRACT 

Wc have developed an efficient photometric technique for identifying young early-type galaxy can- 
didates using a combination of photometric redshifts, spectral-type classification, and optical/near- 
infrared colors. Applying our technique to the GOODS HST/ACS and VLT/ISAAC data we have 
selected a complete and homogeneous sample of young elliptical candidates among early-type field 
galaxies. The distribution of structural parameters for these candidates shows that their selection, 
which is based on early spectral types, is fully consistent with early morphological types. We investi- 
gate the evolution of their luminosities and colors as a function of redshift and galaxy mass and find 
evidence for an increasing starburst mass fraction in these young early-type galaxy candidates at higher 
redshifts, which we interpret in terms of massive field galaxies experiencing more massive/intense star- 
bursts at higher redshifts. Moreover, we find indications for a systematically larger young elliptical 
fraction among sub-L*/2 early- type galaxies compared to their brighter counterparts. The total frac- 
tion among the field early-type galaxies increases with redshift, irrespective of galaxy luminosity. Our 
results are most consistent with galaxy formation scenarios in which stars in massive early-type field 
galaxies are assembled earlier than in their low-mass counterparts. 
Subject headings: galaxies: K+A galaxies — galaxies: formation and evolution 



1. INTRODUCTION 

The last star-formation burst in galaxies defines their 
photometric appearance and is an important diagnostic 
of galaxy formation and evolution. The class of rejuve- 
nated early-type galaxies includes so-called 'K+A' galax- 
ies^ which are objects that show spectroscopic signa- 
tures of old (K-type) and young (A-type) stellar popula- 
tions (Dressier & Gunn 1983). The lack of emission lines 
in their spectra indicates that star-formation processes 
abruptly ceased ~ 1 Gyr ago, followed by a quiescent 
evolution into normal early-type galaxies (e.g. Couch & 
Sharpies 1987). 

Several scenarios have been suggested to explain this 
observed rejuvenation phenomenon: (1) galaxy-galaxy 
mergers (e.g. Zabludoffet al. 1996), (2) interactions of in- 
falling galaxies with the intracluster medium (e.g. Gunn 
& Gott 1972; Bothun & Dressier 1986), (3) galaxy ha- 
rassment (e.g. Moore et al. 1996, 1998), (4) tidally in- 
duced star formation (e.g. Byrd & Valtonen 1990), and 
(5) dusty starbursts (e.g. Poggianti et al. 1999; Pog- 
gianti & Wu 2000). While tidal interaction and galaxy 
harassment are important events in the overall context 
of galaxy formation and evolution (e.g, Bekki et al. 
2005), they do not provide an obvious explanation for 
the sudden halt of star formation without invoking ad- 
ditional processes. Ram-pressure stripping and mergers 
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Historically, these galaxies were named "E-|-As" because of 
their early-type morijhology appearance on the very first images. 
However, subsequent studies have shown that many of these galax- 
ies contain significant disk components, so that the "E" becomes 
unjustified and is substituted by "K", owing to the spectral type 
of the underlying old stellar population (Pranx 1993). 



are currently seen as the most viable scenarios to de- 
scribe the ignition and sudden cessation of star forma- 
tion in young early-type galaxies (e.g. Rose et al. 2001; 
Bekki & Couch 2003; Shioya et al. 2002, 2004; Quintero 
et al. 2004). Moreover, recent radio observations appear 
to exclude a dusty starburst scenario in general, as the 
measured radio fiuxes are inconsistent with high star for- 
mation rates (i.e. -^^ 10 — 100 Mq a^^) for virtually all 
local post-starburst early-type galaxies (Miller & Owen 
2001; Goto 2004). However, Small et al. (1999) does find 
evidence for faint radio fluxes in spectroscopically con- 
firmed post-starburst galaxies in distant dense clusters 
at zKi 0.4. 

K+A fractions are known to vary significantly with 
both redshift and environmental density. Several studies 

have revealed that the brightest K+A's in nearby clus- 
ters are sub-L* systems (Caldwell et al. 1999; Poggianti 
et al. 2004), while at intermediate redshifts (z > 0.8) 
K+As are found with L « 3L* (Dressier et al. 1999; 
Tran et al. 2003), suggesting an increase of K+A clus- 
ter galaxy mass with redshift. So far, the perhaps most 
unbiased comparison between the K+A fractions in field 
and cluster environments was done by Tran et al. (2003, 
2004) who found a field K+A fraction of ~ 3% and a 
factor ~ 4 higher fraction in cluster environment for the 
redshift range 0.3 < z < 1. Tran et al. further noted 
that the field K+A fraction shows strong fluctuations 
and is sensitive to selection criteria. Furthermore, the 
relatively short lifetime of the K+A signature (~ 1 Gyr, 
Couch & Sharpies 1987; Belloni et al. 1995; Barger et 
al. 1996) implies that these galaxies are rare. This indi- 
cates the need for wide-area surveys to identify statis- 
tically signiflcant samples of young early-type galaxies 
and determine their number fractions in both field and 
cluster environments. If their selection is performed con- 
sistently between datasets and galaxy formation models, 
their abundance may be a sensitive diagnostic for star- 
formation activity in field and cluster galaxies as a func- 
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tion of redshift, making these rejuvenated objects useful 
tools to constrain hierarchical galaxy formation scenar- 
ios. 

Keeping in mind that medium-resolution spectroscopy 
of all galaxies in large-area surveys is very expen- 
sive in terms of telescope time, the aim of this pa- 
per is to develop a photometric technique to identify 
young early-type galaxy candidates (hereafter yE can- 
didates) for spectroscopic follow-up using a combination 
of photometric redshifts, spectral-type classification, and 
optical/near-infrared color-color diagrams. Applying this 
technique to the GOODS/CDF-S datasct (Giavalisco ct 
al. 2004) , we select a complete and homogeneous sample 
of yE candidates in the redshift range < z < 1 . Finally, 
we investigate correlations between the colors of yE can- 
didate galaxies with redshift and galaxy mass, which, 
if confirmed by spectroscopy, may put strong constraints 
on galaxy formation models. We use the standard ACDM 
cosmology parameters flu = 0.3, = 0.7, and Hq = 70 
km Mpc""'^ throughout this work. 

2. GOODS/CDF-S DATA 

One of the two fields covered by the Great Observato- 
ries Origins Deep Survey (GOODS) Treasury Program 
is the Chandra Deep Field South (CDF-S) with optical 
HST/ACS and ground-based VLT/ISAAC near-infrared 
photometry (Giavalisco et al. 2004). These data span the 
wavelength range from near-UV to the near-infrared and 
provide ACS imaging with very high spatial resolution 
(0.03") for ~ 86000 galaxies down to V ^ 27.5. Object 
detection and photometry was performed with SExtrac- 
tor (Bertin & Arnouts 1996) after the images in all pass- 
bands were convolved to a common FWHMf«0.45". The 
final magnitudes are the MAG_AUTO SExtractor mag- 
nitudes, which measure the total light. 

Using the multi-waveband data available for galaxies 
in this field, photometric redshifts were measured follow- 
ing the technique presented in Mobasher et al. (2007). 
Briefly, this approach compares the template Spectral 
Energy Distributions (SEDs) for difl'erent spectral types 
of galaxies, shifted in redshift space, with the observed 
SED for galaxies detected in GOODS (see also Sect. 4.1). 

fits are performed and the redshift and spectral type 
correspoiidiug to the minimum value are associated 
to those for the galaxy in question. The template spec- 
tral types consist of elliptical, Sbc, Scd, and Im-type 
galaxies from Coleman et al. (1980) and starburst tem- 
plates from Kinney et al. (1996). Bayesian priors based 
on observed galaxy luminosity functions are used (as 
detailed in Mobasher et al. 2007). Comparison with 
galaxies with available spectroscopic redshifts yields a 
photometric redshift accuracy of a{Az) = 0.03 where 

A^ (^phot ^spcc)/(l ^" ^spcc). 

In order to check the accuracy of our spectral types, we 
compare these with the morphology of galaxies measured 
from the GOODS-S ACS images. The morphologies are 
estimated using the concentration and asymmetry pa- 
rameters (Conselice 1997; Abraham et al. 1996), mea- 
sured for individual galaxies (B. Mobasher- private c;oni- 
munication). We found a very good correlation between 
the morphology parameters and spectral types, particu- 
larly for elliptical galaxies. Therefore, we have reliable 
identification of elliptical galaxies, based on their spec- 
tral types, which constitutes the main selection criterion 



for the present study. 

We use the information on redshifts and spectral types 
to compute if-corrections for all galaxies in our sample. 
Throughout the rest of this work we use AB magnitudes 
(Oke 1974)^, and correct for Galactic foreground redden- 
ing Eb-v = 0.008 derived from the Schlegel et al. (1998) 
maps. We select galaxies brighter than Mx = — 19.0 mag 
from the initial GOODS/CDF-S sample which have a 
mean photometric redshift accuracy of Az < 0.15 cor- 
responding to a 95% confidence interval. This photo- 
metric redshift uncertainty limit is introduced in order 
to keep the average color uncertainty below 0.15 mag 
due to k-corrections (see also Sect. 3). The final sam- 
ple is volume-selected and complete to ^ « 1. We ex- 
clude objects where the photometry was compromised by 
saturated pixels, corrupted or truncated isophotal aper- 
tures, and where the extraction algorithm experienced 
problems during deblending or the photometric errors 
were larger than Am = 0.15 mag in any filter. The fi- 
nal catalog includes objects with HST/ACS BViz and 
VLT/ISAAC JHKs magnitudes'". 

3. STELLAR POPULATION MODELS 
3.1. Diagnostic Colors 

Based on the Bruzual & Chariot (2003) simple stellar 
population (SSP) models and their 2003 version of the 
GAL AXE V stellar population synthesis code, we com- 
pute optical/near- infrared photometric colors for com- 
posite stellar populations (CSP). We use the Salpeter 
IMF with limits at 0.1 and lOOM©. By comparison with 
the control sample of spectroscopically confirmed K-|-As, 
we identify the best diagnostic combination optical/near- 
infrared colors to select yE candidates. 

Given the availability of many filters in our GOODS 
dataset {BVizJHKg) we investigate several filter combi- 
nations that may most efficiently reduce the age-metal- 
licity-extinction degeneracy in the color-color plane. Due 
to the presence of a hotter ^ 1 Gyr young stellar pop- 
ulation, the B-band fiux of yE galaxies is expected to 
be enhanced relative to older (» 1 Gyr) stellar popula- 
tions. For slightly younger stellar populations with ages 
~ 0.4 - 0.8 Gyr the Asymptotic Giant Branch (AGB) 
is densely populated with relatively cool stars (Renzini 
1981) causing a stronger output in the near-infrared K- 
band (Perssou et al. 1983). Since a small population of 
luminous thermally-pulsing AGB stars (TP-AGB) can 
significantly contribute to the integrated light of galaxies, 
the consideration of such short-lived evolutionary phases 
in SSP models is essential (Maraston 1998; Mouhcine & 
Langon 2003). Past this so-called AGB-phase transition, 
near-infrared light is primarily sensitive to the mean tem- 
perature of the red giant branch, which is mainly driven 
by the luminosity-weighted mean metallicity of compos- 
ite stellar populations. 

In Figure 1 we show the effective filter throughputs 
for HST/ACS and the VLT/ISAAC instrument together 
with spectral energy distributions (SEDs) for single-burst 
stellar populations for ages t = 290 Myr, 640 Myr, 1.4 
Gyr, 2.5 Gyr, and 11 Gyr at solar metallicity. For illus- 

® AB magnitudes arc defined by rriAB = —2.5 log /i/— 48.6, where 
/,/ is in ergs s^'^ cui^'^ Hz^'^. Alternatively, one can write tuab = 
— 2.51og/^ -|- 8.9, where is the flux in Jy. 

In the following we shall refer to the Kg filter as the K band. 
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Fig. 1. — Comparison of restframe solar-metallicity SEDs for single-burst stellar populations with ages t = 290 Myr, 640 Myr, 1.4 Gyr, 
2.5 Gyr, and 11 Gyr. A dotted line illustrates the SED of a solar-metallicity stellar population with a constant star formation rate after 
6 Gyr. All SEDs were taken from Bruzual & Chariot (2003). Total filter throughputs for ACS and the ISAAC camera, including mirror 
reflectivity, window transmission, and detector quantum efficiency, are shown at the top of the panel. 



tration purposes all SEDs were normalized to the flux at 
1.3 ^m. It is clear from the figure that the B band con- 
tains the most age-sensitive flux and that the B—J color 
is a very good age indicator. It is instructive to note that 
the Lyman break is a sensitive indicator to identify con- 
stantly star-forming or very young post-starburst galax- 
ies (i.e. t < 0.5 Gyr, see dotted SED in Fig. 1). How- 
ever, the Balmer break is the more sensitive proxy to 
identify post-starburst galaxies with intermediate ages 
(i.e. t > 0.5 Gyr), such as K-l-A galaxies. To find a 
good metallicity indicator we need to move to the near- 
infrared, past the z-band, where the continuum slope 
shows little impact of varying age. Here, the i?-band con- 
tinuum suffers from variations due to strong water vapor 
absorption, which is sensitive to the mean temperature of 
the stellar population (i.e. age-sensitivity). However, the 
relative continuum fluxes in J and K show little variation 
with age, and we deflne J—K as our best metallicity indi- 
cator which exhibits only modest age sensitivity during 
the AGB-phase transition (see also Ferraro et al. 2000). 

The best choice of the most age-sensitive and the 
most metallicity-sensitive color depends also on data 
quality. For reference, the average color uncertainty in- 
troduced by the photometric redshift error is < 0.15 
mag in both B — J and J — K for all galaxy types at 
z < 1. Given the mean uncertainty of all color combi- 
nations of our data, the most age-sensitive color is B — J 
and the most metallicity-sensitive color is J — K, which 
is virtually independent of dust reddening. This BJK 
combination optimizes the selection of post-starburst 
intermediate-age stellar populations at redshifts z < 
1. For higher redshifts, the age-metallicity degeneracy 
dilutes the age/metallicity sensitivity of this color com- 
bination as the restframe flux at A < 1 /im moves into the 
J band, and a different filter combination is required. 



A comparable photometric selection technique, al- 
though specifically targeting star-forming galaxies at 
higher redshifts, was described by Daddi et al. 
(2004). Similar techniques to detect intermediate-age 
globular clusters (i.e. single-burst stellar populations) in 
nearby galaxies have been recently used by Goudfrooij 
et al. (2001), Puzia et al. (2002), and Hempel & Kissler- 
Patig (2004). 

3.2. Spectroscopically confirmed K+As 

To test our photometric selection technique we use 
the K-l-A galaxy sample of Galaz (2000). This data set 
contains spectroscopically confirmed K-|-As in nearby 
(z « 0.05) and distant (z « 0.3) clusters and in the 
nearby field (z « 0.1). All galaxies have been classified 
as K-pAs (Couch & Sharpies 1987; Franx 1993; Caldwell 
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Fig. 2. — Comparison of B magnitudes for spectroscopically 
confirmed K-l-A galaxies taken from Galaz (2000) with those taken 
from the APM catalog. A vertical solid line shows the mean formal 
uncertainty. 
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Fig. 3. — Left panel: Comparison of (B — J)ab vs. (J— ^)ab restframe colors of spectroscopically confirmed K+A galaxies (open 
symbols) with SSP models of Bruzual & Chariot (2003). Open squares mark cluster K-f A galaxies, open circles are field K-(-As. Isochrones 
{solid lines) are plotted for ages t = 200, 500, 700 Myr, and 1, 1.3, 2, 5, and 15 Gyr. Iso-metallicity lines (dotted lines) are shown for 
the metallicities [Z/H] = —2.3, —1.7, —0.7, —0.4, 0.0, and -1-0.4 dex. Right panel: Same data, but this time with SSP models of Anders & 
Fritze-v. Alvensleben (2003). Here isochrones {solid lines) were plotted for the ages t = 200, 500, and 700 Myr, and 1, 1.3, 2, 5, and 14 Gyr. 
Iso-metallicity tracks (dotted lines) are indicated for [Z/H] = -1.7,-0.7,-0.4,0.0, and +0.4 dex. Each panel has the extinction vector for 
Eg_v = 0.2 mag indicated in the upper left corner. Mean uncertainties of the Galaz (2000) data are shown in the upper right corner. 



& Rose 1997; ZabludofT et al. 1996), without any indica- 
tions of emission lines. The near-infrared photometry of 
the Galaz study was performed with the Las Campanas 
NICM0S3 Hg:Cd:Te detector in the CIT filter system 
(Persson et al. 1998). We compute the transformations 
to the (Bessell & Brett 1988, hereafter BB) standard sys- 
tem using Bruzual & Chariot (2003) SEDs with ages from 
0.2 Gyr to 15 Gyr and [Fe/H] values from -2.25 to +0A 
dex, and find the following (BB)-(CIT) transformations 

Jbb = JciT - 0.0045 -I- 0.0102 (J - K)cit 
(J - K)bb = 0.0147 + 0.9835 (J - K)cit 

with RMS scatter of 0.0009 mag and 0.007 mag, respec- 
tively. Galaz' Table 2 lists Johnson B magnitudes as re- 
ported in NED^ at the time, likely coming from various 
original sources. We compare these B magnitudes with 
those in the APM catalog^ (Maddox et al. 1990) and 
show the comparison plot in Figure 2. We find signifi- 
cant offsets for the optical photometry between the Galaz 
and the APM catalog. Since the magnitudes in the APM 
catalog were derived in an internally consistent way, we 
adopt the APM B magnitudes for further analysis in 
this paper. The mean uncertainties in BJK colors of the 
comparison sample are 0.12, 0.10, and 0.18 mag, respec- 
tively. All optical and near-infrared magnitudes from the 
GOODS catalog and the comparison sample represent 
total magnitudes. 

3.3. Confirmed K+As in the BJK color-color plane 

In Figure 3 we compare {B — J)xb vs. {J—K)p^b rest- 
frame colors of spectroscopically confirmed K-f A galax- 

* the NASA/IPAC Extragalactic Database, 
http:/ /nedwww. ipac.caltech.edu 

^ http:/ /www.ast.cam.ac.uk/'apmcat 



ies with SSP model predictions from two different groups: 
Bruzual & Chariot (2003, hereafter BC03) and Anders & 
Fritze-v. Alvensleben (2003, hereafter AF03). The BC03 
predictions were calculated in the AB system with the 
same filter transmission functions as the GOODS/CDF-S 
data. For the AF03 models, which are tabulated in Vega- 
based magnitudes using the BB filter system, we follow 
the Vega-AB transformations: i^AB — 1.891 -f i^vcga, 
Jab = 0.91 -I- Jvega, and Bab = -0.105 -I- i?vcga. As ex- 
pected, the location of confirmed K-|-As in the color-color 
plane corresponds to old to intermediate ages. We find 
that the absolute age calibration differs between the two 
models. The BC03 models predict (luminosity- weighted) 
ages between ^ 15 and ^2 Gyr for most K-|-As, while the 
models of AF03 predict younger ages in the range ^ 1—5 
Gyr. The metallicity scales on the other hand are very 
similar in both models. The majority of confirmed K-f A 
galaxies have luminosity-weighted metallicities between 
[Z/H] —0.7 and 0.0 dex, with few outHers at super- 
solar metallicities, which might be reddened objects. 

3.4. Composite Stellar Populations in the BJK 
Color-Color Plane 

Several groups compared the mean colors of K-f A 
galaxies using optical and near-infrared photometry with 
population synthesis models (e.g. Newberry et al. 1990; 
Belloni et al. 1995; Barger et al. 1996; Shioya et al. 2002; 
Balogh et al. 2005). In this work we combine spectral- 
type fitting with model predictions and investigate the 
infiuence of different metallicities and mass fractions of 
the starburst population on colors and luminosities of 
post-star burst galaxies in the optical/near-infrared color- 
color plane and in color-magnitude diagrams. 

Figure 4 shows the evolution of BC03 isochrones where 
a starburst is ignited after 10 Gyr of passive SSP evolu- 
tion. Assuming a minor merger, we assign 10% of the 
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Fig. 4. — Left panel: Evolution of (B— J)ab vs. {J — K)\^ restframe colors for an SSP with a starburst that ignites after 10 Gyr 
of passive evolution and decays for another 5 Gyr. Two set of isochrones with ages t = 10.1, 10.2, 10.5, 11, 11.5, and 14.9 Gyr are shown 
for the metallicities [Z/H] = 0.4 (solid lines) and —2.3 dex {dot-dashed lines) of the starburst stellar population. The starburst strength 
amounts to 10% of the total final stellar mass. Right panel: Evolution of colors for different starburst strengths / masses for a solar metallicity 
burst. Ischrones for ages t = 10.2, 10.5, 11.5, and 14.9 Gyr are plotted with different line types (solid, dashed, and dotted lines) indicating 
the dependence on starburst mass Mburst = 0.05, 0.1, and 0.5 Afjotal- For reference we overplotted in both panels the SSP model grid from 
the left panel of Figure 3. All predictions were computed using the SSP models of Bruzual & Chariot (2003). Spectroscopically confirmed 
K+A galaxies are as in Figure 3. 



final stellar mass to the starburst population. We adopt 
two metallicities Z = 0.04 (twice solar metallicity) and 
Z = 10""^ ([Z/H] = -2.3 dex) for the starburst stellar 
population, which are superposed on top of the full range 
of metallicities ([Z/H] = -2.3 to -^0.4 dex.) of the un- 
derlying, older stellar population. Note the significant 
difference in B~J colors during the first ~1.5 Gyr after 
the starburst. If determined relative to passively evolv- 
ing galaxies, this color offset is a sensitive post-starburst 
indicator and will be discussed below in more detail. 

In the right panel of Figure 4 we show the influence of 
changing starburst mass fraction on the post-starburst 
galaxy colors. This time we choose solar metallicity for 
the starburst with mass fractions are 5%, 10% (minor 
merger), and 50% (major merger) of the final total stel- 
lar mass. The Figure illustrates that the higher the star- 
burst mass fraction the bluer the B — J color is dur- 
ing the post-starburst phase, in agreement with previ- 
ous studies (e.g. Barger et al. 1996). In addition, there is 
a mass fraction-metallicity degeneracy of post-starburst 
near-infrared colors. Metal-rich major starbursts can ex- 
hibit the same J — K colors as metal-poor minor star- 
bursts. This degeneracy is most prominent for interme- 
diate metallicities in the range —0.7 < [Z/H] < 0. For 
more extreme metallicities the burst strength can be 
constrained with increasing confidence. In other words, 
yE galaxies with starburst mass fractions < 10% and 
J — K > are likely to have experienced a recent 
starburst with super-solar metallicity. The same mass- 
fraction limit and J — K < —0.3 mag indicate sub-solar 
starburst metallicities. 

In summary, B—J is mainly degenerate in age and star- 
burst mass fraction, J—K is mainly degenerate in metal- 
licity and starburst mass fraction. We conclude that the 
B — J vs. J—K color-color plane can be used to identify 



post-starburst galaxies. But only extreme metallicities 
of the post-starburst stellar population can be narrowed 
down robustly. The color-color plane is not suitable to 
derive starburst mass fractions for intermediate metal- 
licities without any other information. However, these 
degeneracies can be partly lifted by including restframe 
K-hand luminosities. Our findings are relatively indepen- 
dent of the choice of SSP models, as the systematics in 
the AF03 models give the same results. 

4. SELECTION OF YE CANDIDATES 

4.1. Spectral Types 

In the following we develop a robust method to iden- 
tify yE galaxy candidates from the GOODS/CDF-S pho- 
tometric catalog using a combination of photometric 
redshifts, spectral-type classification, and optical/near- 
infrared colors. The photometric redshifts and spectral 
types are measured by using template SEDs, shifting 
them in redshift steps and performing fits to the ob- 
served SEDs at each step. The templates corresponding 
to the minimum values are selected, with its spec- 
tral type and redshift assigned to the galaxy in question. 
We use the observed templates for elliptical, Sbc, Scd 
and Im-type galaxies from Coleman et al. (1980) and the 
starburst SEDs from Kinney et al. (1996). In order to 
increase the resolution in spectral types, we also divide 
the interval between adjacent spectral classes into two 
intermediate bins. For instance, between the two dis- 
crete spectral types for ellipticals (T = 1) and Sbc-type 
spectra (T = 2) we include two discrete sub-categories 
(T = l5,l|) that are linear mixtures between the two 
boundary types. This allows a more continuous spectral 
classification across the Hubble types while not causing 
degeneracy. The choice of steps of 1 /3 unit was made as a 
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Fig. 5. — Distribution of spectral types T that were assigned 
by the photometric redshift routine for all objects in our selected 
GOODS/CDF-S sample (for details on the photometric redshift fit- 
ting routine we refer the reader to Mobasher et al. 2004, 2007). The 
locations of individual spectral types are labeled. A sub-population 
of several objects have a slightly earlier spectral type (T = 1 1 ) than 
passively evolving galaxies (T = 1). These are our yE candidates 
and are marked in the plot with a hatched bin. 

compromise between the goal to increase the resolution of 
spectral types on one side and the goal to avoid creating 
degeneracies in the phot-z routine on the other. Details 
of the photometric redshift and spectral-type measure- 
ment are presented in Mobasher et al. (2004, 2007). We 
plot the distribution of the discrete spectral types in Fig- 
ure 5 and define all objects that were assigned slightly 
earlier spectral types {T = 1 ^ ) than passively evolving 
galaxies (T = 1) as yE candidates, i.e. objects which 
are bluer than normal ellipticals but redder than normal 
early-type spirals. Their location in Figure 5 is marked 
by a hatched histogram. Because later galaxy types can 
exhibit similar BJK colors as genuine yE galaxies, we 
use the spectral-type information to exclude late-type 
star-forming galaxies and to include only galaxy candi- 
dates with early spectral types in the following analy- 
sis. The distribution of photometric redshifts of all ob- 
jects is shown in Figure 6. Most of our yE candidates 
have redshifts between z w 0.2 and 0.8, with few objects 
at slightly higher redshifts z k, 1.0. 

4.2. Morphologies 

To investigate the yE candidate morphologies we cross- 
correlate the GOODS/CDF-S sample with the struc- 
tural parameter classification of Ravindranath et al. 
(2006). Based on the BViz HST/ACS images the au- 
thors computed CAS parameters ('asymmetry', 'concen- 
tration', and 'smoothness' defined in Consclicc 2003) for 
all our sample sources. Because of the redshift distribu- 
tion of our sources each ACS filter probes a different 
wavelength and therefore a different morphological (stel- 
lar) component in a galaxy. To allow a comparison of 
structural parameters that is independent of redshift we 
compute the redshift-corrected restframe i3-band mor- 
phologies by linearly interpolating the CAS parameters 
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Fig. 6. — Distribution of photometric redshifts for all selected 
objects in our GOODS/CDF-S photometric dataset (for details on 
the photometric redshift routine see Mobasher et al. 2004). The 
parametrization by spectral-type is identical to the one in Figure 5. 

between the pivot wavelengths of each ACS filter (see 
Sirianni et al. 2005) according to the redshift of each 
source. In Figure 7 we show the distribution of rest- 
frame _B-band CAS parameters for all sample galaxies 
sub-divided by spectral- type as in Figure 5. We find that 
the overall spectral-type classification of our yE candi- 
dates is consistent with early- type morphologies, which 
confirms the robustness of our procedure to select early- 
type post-starburst galaxies. In particular, the distribu- 
tion of the CAS 'concentration' parameter indicates that 
very few yE candidate galaxies with later-type morpholo- 
gies were selected by our photometric selection. 

4.3. BJK Color-Color Diagram 

The left panel of Figure 8 shows the GOODS/CDF-S 
data in the BJK color-color plane. Different symbols 
depict different spectral types as assigned to each ob- 
ject by the photometric redshift fitting routine. Most 
objects with SEDs consistent with passive early-type 
galaxies clump around restframe (_B — J)ab ~ 1-9 and 
(J — iir)AB ~ 0. Objects with shghtly later SED types 
(our yE candidates) scatter towards bluer B—J colors and 
have intermediate J—K colors that place their mean locus 
between spiral and star-forming galaxies. We find that 
the spectral-type selection is consistent with the BJK 
color selection and morphological classification. More- 
over, most yE candidates are consistent with the pre- 
viously computed composite stellar population models 
for post-starburst galaxies. Only one yE candidate at 
{B — J)ab ~ 1-6 and {J—K)ab ~ 0.35 deviates signifi- 
cantly from the model predictions and might be a dusty 
starburst. For reference, we plot the control sample of 
spectroscopically confirmed K-|-As (open symbols). 

To provide a comparison independent of SSP model 
predictions we show in the right panel of Figure 8 
nearby early-type (taken from Michard 2005) and late- 
type galaxies (taken from Perez-Gonzalez et al. 2003a) 
in the same color-color plane as the GOODS/CDF-S 



Evidence for the Evolution of Young Early-Type Galaxies in the GOODS/CDF-S Field 




nj] E/K+A -z 



3 4 
Concentration 




0.2 0.4 
Asymmetry 




0.2 0.4 
Smoothness 



Fig. 7. — Distribution of restframe S-band CAS parameters (concentration, asymmetry, and smoothness, see Conselice 2003) taken from 
Ravindranath et al. (2006) for all selected objects in our GOODS/CDF-S photometric dataset parameterized by their spectral type as in 
Figure 5. Note the strong separation of yE candidates in the concentration parameter [left panel) from galaxies with later spectral types. 
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Comparison between photometric and spectroscopic 
redshifts. 
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data. We find a good agreement between the early-type 
galaxies in GOODS/CDF-S and the local galaxy sam- 
ple. However, the smaller photometric errors of the local 
galaxy sample allow a more clearly detected offset of the 
spectroscopically confirmed K-t-A galaxies towards bluer 
B—J colors (i.e., younger ages) compared to the local 
early- type galaxies, which line up in a tight sequence. 
We also notice that fewer spiral (squares) and Irr galax- 
ies (triangles) lie outside of the SSP grid which is likely 
again, partly due to smaller photometric errors of the 
nearby sample. We investigate whether this disagreement 
may be due to systematics in the photometric redshift 
determination, but find no differences in the error dis- 
tributions between spiral and irregular galaxies, in par- 
ticular m. J — K that appears to be primarily responsi- 
ble for the scatter. A matching between the photometric 
redshifts of our sample and currently available spectro- 
scopic redshifts from the literature (Le Fevre et al. 2004) 
results in good agreement (see Table 1). Although we 
compare restframe colors in both samples, the galaxies 
in our GOODS/CDF-S data are more distant than the 
nearby comparison sample. Evolutionary factors, such 
as possibly enhanced dust fractions, may be responsible 
for the enhanced scatter in our data. Deep mid-IR imag- 
ing and optical to infrared spectroscopy should be able 
to resolve this issue. 

We define the population of yE galaxy candidates 
as those objects that have slightly later spectral types 
(T = li, see Fig. 5) and BJK colors that are inconsis- 
tent with passively evolving (reddest) early-type galax- 
ies. The selected yE candidates are marked in Figure 8 by 



red circles with central black dots. This definition might 
exclude some genuine K-|-As since the control sample 
contains some spectroscopically confirmed K-l-A galaxies 
(in galaxy clusters) with BJK colors consistent with pas- 
sively evolving early- type galaxies. Taking into account 
the average photometric error, this fraction is expected 
to be less than ^ 15%. We point out that even with 
the color, spectral-type, and morphology selection, our 
yE candidate sample may be contaminated by early-type 
spiral galaxies that can mimic the morphologies of ellipti- 
cal galaxies (e.g., when in a largely face-on configuration) 
and the SEDs of yE candidates because of their low, but 
continuous star formation rates. Based on the local com- 
parison sample we estimate this contamination fraction 
between ~ 10 — 30%, depending on the exact color selec- 
tion. Again, deep follow-up spectroscopy should be able 
to address this contamination fraction. 

5. DISCUSSION 

5.1. Color- Magnitude Diagrams 

Post-starburst galaxies tend to be brighter in the opti- 
cal than their passive counterparts (e.g. Tran et al. 2003), 
which is likely due to their additional ~ 1 Gyr old stel- 
lar population. However, optical photometry is not able 
to provide strong constraints on how the stellar masses 
of these rejuvenated galaxies compare with those of pas- 
sively evolving galaxies. The K-hand stellar mass-to-light 
ratio changes by a factor of ~ 3 from '--^ 1.5 Gyr to 13 Gyr 
(see also Drory et al. 2004), which is about an order of 
magnitude smaller a change than in the optical. M/Lk 
changes by about a factor 1.5 between 500 Myr and 1.5 
Gyr (the K-l-A phase). Hence, relative to optical colors a 
Mk selection defines the best approximation to a stellar- 
mass selected dataset. 

In Figure 9 we show color-magnitude (CM) diagrams 
for all early-type galaxies (T < 1.5) and the selected 
yE candidates as well as the K-t-A control sample and 
nearby ellipticals. Overplotted are color-magnitude rela- 
tions for early-type galaxies taken from the large galaxy 
survey of Mobasher et al. (1986), in good agreement 
with the GOODS data and the nearby elliptical galaxy 
sample. We find no significant difference between red- 
sequence galaxies and yE candidates in the near-infrared 
J—K vs. Mk diagram, but we detect a significant offset 
in the B^J vs. Mk diagram where the yE candidates 
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Fig. 8. — Left panel: BJK color-color diagram (using de-reddened rest-frame AB magnitudes) of galaxies in the GOODS/CDF-S Field. 
The data symbols are parameterized by the galaxy spectral types, and are early-type {filled circles), spiral {filled squares), irregular {filled 
triangles), and starbursts {stars). yE galaxy candidates are marked by small dark dots. Spectroscopically confirmed K+A galaxies taken 
from Galaz (2000) are indicated by open symbols, where open squares mark cluster K-l-A galaxies and open circles are field K-I-As. Average 
total photometric errors, including uncertainties of the photometry and the photometric redshift determination, are indicated in the upper 
right corner for yE candidates {solid cross) and the mean photometric uncertainties for confirmed K-l-A galaxies {dotted cross). Right 
panel: BJK color-color diagram for local galaxies. Open circles are nearby ellipticals from Michard (2005), open squares and triangles are 
Sab and Scd spirals, while stars show local starburst galaxies, both taken from Perez-Gonzalez et al. (2003a). The same spectroscopically 
confirmed K-l-A sample as in the left panel is overplotted for comparison. The solid error bar shows average uncertainties for the elliptical 
sample, while the dotted cross is for the later-type galaxies. SSP models in both panels are identical to those in the left panel of Figure 4. 



have significantly bluer colors at the same Mk (stellar 
mass) compared to the red sequence. 

We show the influence of changing luminosity- weighted 
age {At), metallicity (AZ), total galaxy mass (AMgai.), 
and starburst-mass fraction (AM*) as vectors in both 
CM diagrams. At and AMgai., and AZ and AM* are 
highly degenerate in the J—K vs. Mk diagram. In fact, 
the At and AMgai. vectors are almost parallel to the red 
sequence (see upper panel in Fig. 9) and have, therefore, 
negligible influence on the scatter in the J—K color, 
which is mainly driven by metallicity and the starburst- 
mass fraction. The relatively small scatter in this dia- 
gram hints at a small scatter in metallicity and starburst- 
mass fraction and/or is a signature of a correlation be- 
tween these parameters. This degeneracy can be partly 
lifted with the B — J vs. Mk diagram (see below). Four 
yE candidates fall significantly off the red sequence in 
the J—K vs. Mk diagram. They are likely to have expe- 
rienced significantly different starburst events than the 
rest of the sample. One possible explanation is that the 
metallicity of their starburst population is significantly 
lower. 

The most influential parameters in the B — J vs. Mk 
diagram are age and the starburst mass fraction (besides 
the total galaxy mass). Both vectors are highly inclined 
relative to the red sequence, which is a graphic explana- 
tion of the scatter in this plot. Metallicity plays a negli- 
gible role. We use the B — J vs. Mk diagram to discuss 
the redshift evolution of our yE candidates in the next 
section. 



5.2. Starburst Mass Evolution 

In the hierarchical merging scenario of galaxy forma- 
tion, massive structures are expected to form on ex- 
tended timescales from smaller sub-units with a consid- 
erable fraction of stars forming relatively recently (e.g. 
Springel et al. 2005). If, in contrast, more massive struc- 
tures form first, as seen in the monolithic collapse sce- 
nario, the majority of stars form at high redshifts {z > 2) 
in intense starbursts (e.g. Larson 1975). The fraction of 
yE signatures among massive early-type galaxies is a 
measure of dissipative merging activity (recently often 
named "wet merging"), which is expected to be a func- 
tion of redshift in the hierarchical merging picture. In 
contrast, no such evolution is expected in the early mono- 
lithic collapse scenario out to z ^ 2. This difference is 
expected to depend on environment and to be more pro- 
nounced in the field (Benson et al. 2002), which we probe 
here with the GOODS dataset. The fraction of yE galax- 
ies and the strength of the yE signatures among massive 
early-type field galaxies as a function of redshift therefore 
puts strong constraints on galaxy formation models. 

To investigate the redshift evolution of our yE candi- 
dates, we subdivide the sample into four redshift bins 
between z = and 1. For each galaxy we determine the 
B—J offset with respect to the red sequence and plot 
the residuals versus their absolute K-hand luminosity in 
Figure 10. Our analysis is limited in redshift space by the 
completeness of our GOODS data in the K band. A lOcr 
point source detection in the VLT /ISAAC data is feasible 
down to TO/f(AB) = 25.1 mag (Giavalisco et al. 2004), 
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Fig. 9. — Optical/near-infrared color-magnitude diagrams of 
early-type galaxies {filled circles) and yE candidates {dotted cir- 
cles). Large open (magenta) symbols mark spectroscopically con- 
firmed K-|-As (Galaz 2000) in the field (circles) and cluster envi- 
ronment {squares). Small open (red) circles are nearby elliptical 
galaxies from Michard (2005), as in the right panel of Figure 8. 
The nearby data were corrected for passive evolution, according 
to Bruzual & Chariot (2003) models, to the reference redshift 
z = 0.5 of our GOODS data, which includes only a correction of 
A{B— J) = 0.2 mag. Arrows indicate the relative changes in At be- 
tween 500 Myr and 1.5 Gyr (short arrow), and 500 Myr and 13 Gyr 
(long arrow), AZ between 0.2Z and solar metallicity, AMg^i. be- 
tween lO^" and 10^1 M0, and AM* between 5% and 50% of the 
starburst mass fraction. Color-magnitude relations for early-type 
galaxies from the survey of Mobasher et al. (1986) are overplot- 
ted as dashed lines. Average photometric errors of yE candidates 
{solid cross) and spectroscopically confirmed K-|-As {dotted cross) 
are plotted in the upper right corner of each panel. 

which translates into Mk, ab = —19.0 mag at a redshift 
of 2; = 1 (indicated by a vertical line in Figure 10). Hence, 
completeness is not an issue for the following analysis. 

We observe that at any given Mk the B—J residuals are 
systematically larger for higher redshift galaxies. Spec- 
troscopically confirmed K-t-As show a similar trend, but 
their redshift range is significantly smaller {z < 0.12) 
than that of our sample {z < 1, see Fig. 6). Note that 
confirmed field K-t-As show larger 6{B—J) than confirmed 
cluster K-f-As, which he closer to the red sequence. We 
also find tentative evidence for a systematic difference 
in {5{J — K)) between confirmed field K+As (open cir- 
cles) and the combination of our yE candidates and con- 
firmed c/wsier K-|-As. It appears that a metallicity differ- 
ence is the most obvious way to explain this offset, given 
that in the more age-sensitive B — J panel the yE candi- 
dates lie among, or are bluer than, the nearby confirmed 
field and cluster K-f As. However, the photometric errors 
of the spectroscopically confirmed K-l-A sample require 
that this intriguing result be confirmed with deep spec- 
troscopy. 

Parameters likely responsible for the evolution in 5{B— 
J) with redshift are (i) a younger age of the post- 
starburst population at higher redshift, or (ii) a larger 
mass fraction of the post-starburst population at higher 




Fig. 10. — Residual B—J colors with respect to the red sequence 
(see Fig. 9) as a function of absolute K band magnitude for early- 
type galaxies. Dark dots mark yE candidates. Open symbols in- 
dicate spectroscopically confirmed K+A as in Fig. 9; taken from 
Galaz (2000). The symbols size is parametrized by redshift. All 
confirmed K-l-A galaxies are located at redshifts z < 0.12. Vec- 
tors indicate the influence of varying age, metallicity, and star- 
burst mass fraction as in Fig. 9. Average photometric errors are 
indicated in the upper left corner of the panel, where the dotted 
cross indicates typical errors of the confirmed K-I-As and the solid 
cross shows errors for yE candidates. A long vertical line shows 
the lOo" points source detection limit at 2 = 1. Small vertical lines 
indicate L* at redshifts 0.2 and 1.0 for field galaxies (Drory et al. 
2003). Note the paucity of yE candidates at magnitudes fainter 
than Mk, AB ~ -21.2 and S{B- J) ab < -0.8 mag. 

redshifts (we established in Section 3 that metallicity has 
a smaller effect on 5{B — J)). First we consider age ef- 
fects. Under the assumption that the yE candidates in- 
deed host 0.5 — 1.5 Gyr old post-starburst populations 
with strong Balmer absorption-line spectra, the vectors 
in Figure 10 show that the difference in d{B—J) between 
yE candidates with z < 0.4 and those with z > 0.6 is 
in principle similar to the effect of age fading from 0.5 
to 1.5 Gyr. However, it is hard to imagine that post- 
starburst populations in yE galaxies at higher redshift 
(be it induced by galaxy interactions or by ram-pressure 
stripping) would be systematically younger than in low- 
redshift ones, given the short duration of the yE phe- 
nomenon relative to the difference in look-back times 
within this redshift range. In the adopted WMAP cos- 
mology, the difference in look-back times between z — 0.2 
and 1.0 is 7.8 — 2.4 = 5.4 Gyr, significantly longer than 
the duration of the yE phenomenon. Instead, one would 
expect the ages of such post-starburst populations to be 
randomly distributed (between 0.5 and 1.5 Gyr) at any 
redshift. We, therefore, suggest that it is more likely that 
the increase of yE candidates' 6{B — J) with redshift is 
primarily due to an increasing mass fraction of the post- 
starburst population. If indeed the yE phenomenon iden- 
tifies an important era in the assembly history of early- 
type galaxies, the lack of low-luminosity {Mfc > —21) 
yE candidates dX z > 0.6 in our GOODS sample would 
constitute strong evidence in favor of galaxy formation 
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scenarios in which more massive early-type galaxies in 
the field are assembled earlier than their low-mass coun- 
terparts, which is in line with the "downsizing" picture 
(e.g. Cowie et al. 1996). 

On the other hand, we currently cannot exclude 
the possibility that the higher-redshift yE candidates 
do not actually host a (one) post-starburst popula- 
tion, but instead have systematically younger overall 
(luminosity-weighted) ages than the lower-redshift yE 
candidates. More frequent bursts at higher redshifts, due 
to more frequent mergers, might be a possible expla- 
nation, which would lead to systematically younger ob- 
served ages. Given the relevance of this issue in terms 
of its power to constrain early-type galaxy formation 
scenarios, we suggest a spectroscopic followup to test 
whether our yE candidates really do contain K-|-A fea- 
tures in their spectra. 

5.3. Evolution of the yE Fraction 

The fraction of post-starburst galaxies in nearby clus- 
ters was reported to increase with decreasing galaxy 
luminosity (e.g. Poggianti et al. 2001; Small et al. 
2001). The morphology bias prevents a direct compar- 
ison between the yE fractions in the cluster and field 
environment. We, therefore, select only early-type galax- 
ies from our sample and investigate the yE fraction in the 
field. A similar trend of an increasing post-starburst frac- 
tion with decreasing galaxy luminosity is present in our 
GOODS sample of field galaxies. We quantify this trend 
by splitting the data at L*/2 along the pure luminosity 
evolution vector and counting the number of passive Es 
and yE candidates. The corresponding values for the di- 
vision for each redshift bin into super- iy*/2 and sub-iy*/2 
systems are = -22.27,-22.43,-22.64, and -22.75 
mag for z = 0.3,0.5,0.7, and 0.9, respectively. Note that 
the samples are complete down to Mb = —16.3 mag and 
Mk = —19.0 mag at z = l, which leaves our sample sen- 
sitive to objects with restframe colors B — K < 2.7 mag, 
but excludes extremely reddened objects, such as EROs 
(e.g. Daddi et al. 2000). However, these galaxies are most 
frequent beyond redshift unity (e.g. Georgakakis et al. 
2006; Simpson et al. 2006) and are expected not to be a 
significant part of our initial sample. To improve sample 
statistics we merge the number counts of the two upper 
and two lower redshift bins, which results in a high and 
low-redshift sample divided at 2; = 0.6. Galaxy number 
counts are given in Table 2. 

We confirm the trend of a higher fraction of post- 
starbursts among early-type galaxies toward less lumi- 
nous objects. In the high-redshift sample we find a yE 
candidate fraction of 86% among all sub-Z/*/2 early-type 
galaxies, whereas only 22% of the more luminous (mas- 
sive) objects are yE c;aiididates. The yE fractions shrink 
towards lower redshifts to 51%. for sub-L*/2 systems and 
11% for their brighter counterparts. We suggest that 
these high yE fractions be confirmed (or refuted) using 
follow-up spectroscopy. 

6. CONCLUSIONS 

Based on the combination of photometric red- 
shifts, spectral- type classification, and optical/near- 
infrared colors, we have selected a sample of young 
early-type galaxy candidates (yE candidates) from 
the GOODS/CDF-S dataset. Our technique relies on 



TABLE 2 

Number counts of E galaxies and yE candidates. 
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> L*/2 






redshift 




Ne NyE 


/<L*/2°' 




z > 0.6 


2 12 


7 2 


0.86 


0.22 


z < 0.6 


18 19 


16 2 


0.51 


0.11 



^Fraction /<£,* /2 of yE galaxy candidates among all sub-L*/2 carly- 
typc galaxies. 

^Fraction />i,* /2 of yE galaxy candidates among all > L*/2 early- 
type galaxies. 



spectral- type fitting and color-color/color- magnitude di- 
agnostic diagrams and provides an efiicient method for 
selecting yE candidates using only photometric infor- 
mation. An analysis of CAS parameters shows that the 
selected yE candidates have early morphological types, 
which confirms that spectral-type and morphological- 
type selection are fully consistent with each other. The 
comparison of population synthesis models and observed 
properties of a sample of spectroscopically confirmed 
K-l-A galaxies provides strong circumstantial evidence 
that the selected candidates axe genuine field early-type 
post-starburst galaxies. 

We study the systematics that drive colors and mag- 
nitudes of post-starburst galaxies in diagnostic diagrams 
using current population synthesis models. Our analy- 
sis reveals evidence for a changing starburst mass frac- 
tion with increasing redshift in the sense that more mas- 
sive/intense starbursts may be responsible for the yE 
signatures of massive field early-type galaxies at higher 
redshifts. Furthermore, we find a higher yE candidate 
fraction in sub-L*/2 early- type galaxies compared to 
their more luminous counterparts. Within the redshift 
range of our data (^ ^ 1) we may be witnessing evi- 
dence for enhanced merging in the field towards higher 
redshifts, especially for high- luminosity galaxies. Similar 
results are obtained by studies of luminosity-weighted 
mean ages of nearby early-type galaxies which indi- 
cate that low-L galaxies show a much greater scatter 
to younger luminosity-weighted mean ages than high-X 
galaxies (e.g. Caldwell et al. 2003; Perez-Gonzalez et al. 
2003b; Sanchez-Blazquez et al. 2006). If follow-up spec- 
troscopy of the yE candidates identified in this paper 
confirms their K-l-A nature, this suggests that the assem- 
bly of stars in high-L early-type field galaxies occurred 
earlier than in their lower-L counterparts. 
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